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Abstract 

The next-to-minimal supersymmetric model (NMSSM) with a large A (the mixing parameter 
between the singlet and doublet Higgs fields) is well motivated since it can significantly push up 
the upper bound on the SM-like Higgs boson mass to solve the little hierarchy problem. In this 
work we examine the current experimental constraints on the NMSSM with a large A, which include 
the direct search for Higgs boson and sparticles at colliders, the indirect constraints from precision 
electroweak measurements, the cosmic dark matter relic density, the muon anomalous magnetic 
moment, as well as the stability of the Higgs potential. We find that, with the increase of A, 
parameters like tan/3. Ma, and M2 are becoming more stringently constrained. It turns out 
that the maximal reach of A is limited by the muon anomalous magnetic moment, and for smuon 
masses of 200 GeV (500 GeV) the parameter space with A ^ 1.5(0.6) is excluded. 
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I. INTRODUCTION 

Since the minimal supersymmetric standard model (MSSM) [l| suffers from the /i-problem 
and the little hierarchy problem, some non-minimal supersymmetric models have recently 
attracted much attention, among which the most intensively studied is the next-to-minimal 
supersymmetric standard model (NMSSM) |^. In the NMSSM there is no dimensionful 
parameters in the supersymmetry-conserving sector and the fi term is dynamically generated 
through the coupling between the two Higgs doublets and a newly introduced singlet Higgs 
field which develops a vacuum expectation value of the order of the SUSY breaking scale. 
The NMSSM provides two ways to alleviate the little hierarchy problem. One is to relax the 
LEP II lower bound on the mass of the SM-like Higgs boson, h, by diluting ZZh coupling 
through the singlet component of h and/or by suppressing the visible decay h —>■ bb through 
introducing new decay of h [4]- The other is to push up the Higgs boson mass with a large 
A, which can be seen from the tree level upper bound of the Higgs boson mass 5| 

mlma. - ^1 cos^ 2(3 + W sin^ 2(3 (1) 

where tan (3 = (Hy) / (Hy) , = (Hu)^ + (H^)'^ and A is the mixing parameter between the 
singlet and doublet Higgs fields defined in Eq.([2]). 

Note that the choice of a large A to solve the little hierarchy may be limited by the 
perturbativity of the theory at the scale A since the value of A is increasing with the energy 
scalelGl. If this scale A is the grand unification (GUT) scale, A should be less than about 0.7 

n 

at weak scale, leading to an upper bound on the Higgs boson mass of about 150 GeV [5|. 
However, the bound on A from the perturbativity consideration can be relaxed by embedding 
the NMSSM in some more complex frameworks. For example, in the Fat Higgs model [tI, by 
completing the NMSSM (or NMSSM-like models) with an appropriate strong dynamics at 
an intermediate scale (much lower than the GUT scale), A can be as large as 2 at weak scale 
and the Higgs boson mass can be pushed up to about 400 GeV. In this work, regardless the 
detailed forms of the ultraviolet physics, we treat the NMSSM as an effective theory and 
examine the current experimental constraints on its parameter space. 

Such phenomenological studies on the Higgs boson and supersymmetry are pressing since 
the mystery of the Higgs sector will be unveiled at the LHC in the near future. If the SM- 
like Higgs boson is discovered with a mass above the MSSM upper bound, the NMSSM (or 



other NMSSM-like models) with a large A, generally called ASUSY 



8|, will be immediately 



favored since it not only inherits all the advantages of the MSSM, such as unifying gauge 
couplings and providing a dark matter candidate, but also is free from the /x-problem and 
the little hierarchy problem. For the phenomenological studies of these models, a primary 
work is to examine the current experimental constraints on their parameter space. 

We note that various constraints on the NMSSM have been studied in the literature, but 
different constraints were considered in different papers. For example, in the authors 
mainly considered the LEP II constraints and put emphasize on small A case. The package 
NMSSMtools encoded various constraints (like the LEP II searches for the Higgs boson, 
the cosmic dark matter and the stability of the Higgs potential), but it is still not complete 
since it does not include the indirect constraints from precision electroweak measurements 
and the muon anomalous magnetic moment. In this work we consider all these constraints 
and especially focus on the case with a large A. As will be shown from our study, with the 
increase of A, the parameter space is getting more stringently constrained. To figure out the 
allowed parameter space is helpful for exploring such low energy supersymmetry at the LHC 
and also may shed some light on constructing the ultraviolet physics from the bottom-up 
view. 

This paper is organized as follows. In Sec. II we briefly describe the structure of the 
NMSSM with emphasis on its difference from the MSSM. In Sec. Ill we summarize the 
constraints considered in this work and briefly discuss their characters. In Sec. IV we scan 
over the NMSSM parameter space and display the region allowed by all these constraints. 
In Sec. V we give our conclusions. 

II. ABOUT THE NMSSM 

The NMSSM extends the matter fields of the MSSM by adding one gauge singlet super- 
field S, and its superpotential takes the form [3j 

W = XsiMHiS + ^kS^ + Y^e,,Q'UHi - Y^e^.Q^M', - Y^SijUEW^ (2) 
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where Q, U and D are squark superfields, L and E are slepton superfields, Hu and Hd are 
Higgs doublet superfields. The soft SUSY breaking terms are given by 

1 



^soft = 7^M2X''X'' + -M^X'X' + mj\Hd\' + mi\H^\' + mi,\S\ 

,2 |/S|2 I _2 irn2 I _2|7^|2 I _2|?|2 , _2 | fi|2 

,1 



+m^|Qp + ml\U\^ + ml\D\^ + ml\L\^ + m||E| 



+{XAxei,HlH',S + h.c.) - {^A^S^ + h.c.) 
+iYuAueijQ'UHi - YdADSijQ'DHi - Y.AESijVEHi + h.c). (3) 

Note that just like the MSSM, the NMSSM has the feature that SUSY breaking induces 
the electroweak symmetry breaking. Before SUSY breaking (i.e. without the soft break- 
ing terms), the Higgs scalars have zero vevs in the supersymmetric vacuum of the scalar 
potential and thus the electroweak symmetry is not broken. After SUSY breaking (i.e. 
with the soft breaking terms), the Higgs scalars develop non-zero vevs in the physical (non- 
super symmetric) vacuum of the scalar potential and hence the electroweak symmetry is 
spontaneously broken and the /i parameter is generated = X{S). Since both the elec- 
troweak symmetry breaking and the /z parameter generation are induced by SUSY breaking, 
their scales should be naturally at the SUSY breaking scale (the scale of soft breaking mass 
parameters). 

The differences of the NMSSM and MSSM come from the Higgs sector and the neutralino 
sector [Si]. In the Higgs sector of the NMSSM there are three CP-even and two CP-odd Higgs 
bosons. In the basis [i?e(if°), -Re(if°), -Re(S')], the mass-squared matrix elements for CP- 
even Higgs bosons are 

-^1,11 = rn\cos^ 13 + m\ srn^ [3 , (4) 

A^l 22 — sin^ (3 + in^z "^"^s^ (5) 

-M|33 = ^misin2 2/3-^t;2sin2/3+^/i(4/€V-AA.), (6) 

M%^2 = (2AV -m| -m^)sin/?cos/?, (7) 

-^1,13 — 2A/xf sin/? — — m^sin2/3cos/3 — K/it>cos/3, (8) 

-^1,23 — 2A/iT;cos/3 — ^m^sin/3sin2/3 — K/iT;sin/3. (9) 
In the basis [A, Jm(S')] with A = cos/? Im{Hl) + sin/? Im{H^), the mass-squared matrix 
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elements for the CP-odd Higgs bosons are 

2/i XAx + Kfi 



M 



2 

PAl 



m 



Ay 



P,22 



PA2 



sin 2/5 A 

3 X^v^ 3 

-A/tf ^ sin 2/9 H ^"^4 sin^ 2(3 + — /iA^, 

2 Afj,^ X 

sin 2/5 — 3Kfiv. 



(10) 

(11) 
(12) 



As shown in Eq. fllOp . we can choose instead of Ax as a free parameter. So compared with 
the MSSM, the NMSSM has three additional parameters: A, k and A^. Conventionally, A is 
chosen to be positive while n and can be either positive or negative. Note that Eqs.dHJ 
and (fT2|) indicate that the parameters A and nfi affect the mixings between doublet and 
singlet Higgs fields, while only affects the squared-mass of the singlet Higgs field. 

In the neutralino sector, the NMSSM predicts one extra neutralino. In the basis 
{—iXi, —1X2, V'u, V's) the neutralino mass matrix is given by 

f Ml sin 6'vi/ sin /5 — sin 6*^/ cos/3 ^ 

M2 —mz cos 9w sin /3 rriz cos 9\y cos P 

-fi -Aw cos /3 . (13) 

— Afsin/? 

V 2f/. ; 

This mass matrix is independent of A^, and the role of A is to introduce the mixings of tps 
with ijj'l and ip^, and kfi is to affect the mass of ips- From Eq. (I9IT2IT3!) one can learn that 
in the limit A, k ^ 0, the singlet field have no mixing with the doublet field and thus is 
decoupled. In this case, the NMSSM can recover the MSSM. 



III. CONSTRAINTS ON THE NMSSM PARAMETERS 



Before we proceed to discuss experimental constraints on the parameters of the NMSSM, 
we take a look at the bounds on A and n from the requirement that the theory should keep 
perturbative under a certain scale A. The renormalization group equations (RGEs) for A 



and K, under the scale A take the following form 
dX X 



dln^ 

dn 
rfln/i 



6k, 
16^ 



(4A^ + 2k^ + 3r/ + 3Y^^ + If - 3^2 - g'^) 
(A^ 



(14) 
(15) 



where g and g' are the SU{2)l and f/(l)y gauge couphngs. These RGEs indicate that the 
values of A and k increase with the energy scale. The requirement of perturbativity till 
the cut-off scale A, i.e., A(A) ^ 27r and k(A) ^ 2-7?, will set upper bounds on A and k at 
weak scale (throughout this paper, without specification all input parameters are defined at 
weak scale). For example, if we assume that new dynamics appears at A = lOTeV, we get 
A^ + ^ 4.2 and for A > 1.5, k must be less than 1.2; while if A is chosen to be the GUT 
scale, a stringent bound A^ + ^ 0.5 is obtained 6|. In our following numerical study we 
let A and k to vary below 2 and 1, respectively, and this corresponds to set A ~ 10 TeV. 
In our study we consider the following constraints on the parameters of the NMSSM: 

(1) Constraints on the neutralino and chargino sector, which include: the bound from 
invisible Z decay V{Z XiXi) < 1-76 MeV; the upper bounds on neutralino pair 
productions at LEP II cr(e+e^ XiX^) < 10^^ pb {i > 1) and cr(e+e^ x'iXj) < 
10^^ pb; and the LEP II bound on the lightest chargino mass m+ > 103.5 GeV. 
These bounds will mainly constrain the parameters Mi, M2 and /i. 



(2) Lower bounds on sparticle masses from LEP and Tevatron experiments 
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me > 73 GeV, > 94 GeV, ruf > 81.9GeV, rriq > 250 GeV, 
mi > 89 GeV, > 95.7 GeV, m-g > 195 GeV, 

where mq denotes the masses for the first two generation squarks. These constraints 
will put lower bounds on the soft breaking masses for sleptons and squarks. 

(3) The LEP II lower bound on the charged Higgs boson mass, mH+ > 78.6 GeV, which 
gives a lower bound on mA through the relation = m\ + M^r — |A^w^. 



(4) Constraints from the direct search for Higgs boson at LEP II 13|], which include 
various channels of Higgs boson productions [l^ . They will constrain the parameters 
m^, tan/5, A as well as the masses and the chiral mixing of top squarks in a complex 
way. 

(5) Constraint from the relic density of cosmic dark matter, i.e. 0.0945 < VlK^ < 0.1287 

n 

|14J |. assuming the lightest neutralino is the dark matter particle. The relic density 



will constrain the parameters Mi, M2, /i, m^, tan/? and A in a complex way 
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(6) Constraint from the stabihty of the Higgs potential, which requires that the physical 
vacuum of the Higgs potential with non- vanishing vevs of Higgs scalars should be lower 
than any local minima. Also, the scale of the Higgs soft breaking parameters should 
not be much higher than the electroweak scale to avoid the fine-tuning problem. Here 
we set 1 TeV as the upper bound of the soft breaking parameters in the Higgs sector. 
This will constrain the parameters tua, /i, ^k, A and tan/5. 

(7) Constraints from precision electroweak observables such as piept, sin^ ^gj^* and M^, 
or their combinations ej(i = 1,2,3) [16]. We require the predicted in the NMSSM 
to be compatible with the LEP/SLD data at 95.6% confidence level or equivalently 



X^/dof < 8.1/3. We take the correlation coefficient of from 17| in calculating x^- 
This requirement will constrain the parameters tan/3, as well as the soft breaking 
parameters in the third generation squark sector. 

(8) Constraint from Ri, = T{Z hh)/T{Z hadrons), whose measured value is Rl^^ = 
0.21629 ± 0.00066 and the SM prediction is = 0.21578 for = 173 GeV [12;]. In 
our analysis we require i?f is within the 2a range of its experimental value. It has 
been shown that the SUSY contribution to Rb might be sizeable for large tan/? jisl ]. 

(9) Constraint from the muon anomalous magnetic momentum a^. Now both the theo- 
retical prediction and the experimental measurement of have reached a remarkable 
precision, but they show a significant deviation a^^^ — aj^^^ = (29.5 ± 8.8) x 10^^° jigj ]. 
In our analysis we require the SUSY effects to account for such deviation at 2a level. 
The character of the SUSY contribution to is that it is suppressed by smuon masses 
but enhanced by tan (3. 



Among the above constraints, (1-6) and (9) have been encoded in the package NMSSM- 



Tools 



10|. In our calculations we extend it by including the constraints (7) and (8) 



The analyfc express,™, of , and „> the NMSSM were given ,„ our reeent work |18|. 
In [18[ we also calculated the NMSSM contribution to (when we started that work, the 
results in |2 



24i\ had not yet published), where we extended the neutralino- and chargino- 
mediated MSSM contributions [2^ to the NMSSM and also considered the contributions 



from the Higgs-mediated diagrams |21j and from the Barr-Zee diagrams 



that our results in [18| agree with those in [23j. 



22] . We checked 
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Note that in our analysis we did not include the constraints from various S-decays 25 1 
because they are dependent on squark flavor mixings and thus involve additional parameters. 



IV. ALLOWED REGIONS OF THE NMSSM PARAMETERS 

In this section, we scan over the NMSSM parameter space to look for the region allowed by 
the constraints in the preceding section. Since we are interested in the parameters sensitive 
to the constraints, we make some assumptions (as conservative as possible) for the other 
parameters such as soft breaking parameters in squark, slepton and gaugino sectors. 

For the parameters in squark sector, we assume the so-called m^"'^ scenario, which can 



maximize the lightest Higgs boson mass 26||. This scenario assumes all the soft breaking 



masses in the squark sector to be degenerate 

M,- = Mq^ = Mu^ = Mn^ (16) 

with i being the generation index. It also assumes the trilinear couplings to be degenerate 
= Ad^ with (A„. — fi cot P)/Mq = 2. We fix Mg = ITeV in our analysis since large Mg 
can not only enhance the lightest Higgs boson mass, but also decrease the contribution of 
the third generation squarks to the electroweak parameters, which has the same sign with 
the Higgs contributions Isj]. For the parameters in slepton sector, we note that the slepton 
masses affect little on the constraints except the muon anomalous magnetic momentum. 
In our calculation we assume all the soft breaking parameters in the slepton sector are 
degenerate and take a value of 200 GeV (we will discuss the effects of its variation). For the 
gaugino mass parameters, we assume the grand unification relation Mi = ^{g''^/g^)M2. 

With the above assumptions, the free parameters are reduced to seven (A, k, A^, tan/3, 
TTiA, /i, M2) and within the capability of our computer to perform a scan. During our scan, 
we first divide the varying range of A into bins with each bin width being 0.1 and then we 
vary the values of other parameters in the following ranges 

-1 < K < 1, l<tan/3<60, 

-1 TeV <A^<1 TeV, 50 GeV < Ma, /x, M2 < 1 TeV. (17) 

With two hundred million samples in each bin and keeping the points satisfying the 
constraints, we finally get the allowed regions of these parameters. Our scan results indicate 
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FIG. 1: The scatter plots are the NMSSM parameters satisfying all the constraints (1-9). The 
curve is the upper bound on tan /? without considering the muon g-2 constraints. 



that the number of the survived samples for A < 0.5 is much larger than that for A > 0.5, 
which means that the parameters for small A are much less constrained than the case with 
large A. Since we are interested in large A, here we only show our scan results for A > 0.5. 

In FigHlwe display the parameters (scatter plots) satisfying all the constraints (1-9) in 
the plane of A versus tan/5. Also, we present a curve which is the upper bound on tan/3 
without considering the muon g-2 constraints. To get this curve, we fix A and scan over 
the parameters in Eg. (1171) . We adopt the important sampling method 27| to optimize the 
varying range of tan /?. 

FigH] shows that the upper bound on tan /? gets stronger as A gets large, and when all the 
constraints are considered, A is upper bounded by about 1.5. The underlying reason for this 
is that the constraints (1-8), especially the constraint (7), have limited the maximal value 
of tan/5, which decreases with the increase of A. Since a large tan/? is needed to explain 
the deviation of the muon g-2, A must terminate at a certain value where the corresponding 
tan/? value is too small to explain the muon g-2. We have checked that the maximal value 
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FIG. 2: Scatter plots of the NMSSM parameters satisfying all the constraints (1-9), displayed in 
different planes. 

of A is dependent on slepton mass. For example, for slepton mass of 100 GeV, 280 GeV and 
500 GeV, the bounds on A are A ^ 2, A ^ 1 and A ^ 0.6, respectively. 

In Figj2] we display the NMSSM parameters satisfying all the constraints in different 
planes. We see that for a large A the parameters m^, /i, M2 and are also bounded 
in a certain region. For A = 1, these bounded regions are 400 GeV ^ Ma ^ 800 GeV, 
150 GeV < /i < 250 GeV, 150 GeV < M2 < 300 GeV and < 600 GeV. 

From the figure of Ma versus A in Fig 12] one can see that the lower bound of Ma increases 
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as A becomes large. The reason is that the LEP II direct search for Higgs boson mainly 
limits the mass and the couplings of the light CP-even Higgs boson whose component is 
dominated by the doublet Higgs field Hu or H^. For tan/5 > 1, this Higgs boson should be 
dominantly composed by if„ field since A^ln is smaller than A^|22' ^^"^ mass is to be 
reduced by the off-diagonal elements 12 and 13- As A gets larger, these off-diagonal 
elements get larger and hence reduce the mass of the light CP-even Higgs boson, which then 
requires a larger to compensate in order to satisfy the LEP II lower bound. 

The figure of /i versus A in Figl2] indicates that with the increase of A, the upper bound 
of n decreases. This is because in the off-diagonal elements A^|i3 and (which reduce 

the light CP-even Higgs boson mass), A is always associated with /x, and to meet the LEP 
II bound a large A must be accompanied by a small /i. 

The figure of M2 versus A in Figl2] shows that M2 is also bounded in a narrow region. 
This is because the relic density of the dark matter correlates the parameters m^, /i, M2, A 
and tan /? in a complex way, and a large value for any of these parameters will limit severely 
the region of other parameters. 

The figure of versus A in FigE] shows that the trilinear soft breaking parameter 
for the singlet field is also limited. This can be understood from the expressions of A^lsa 
and A^p22- The stability of the Higgs potential requires both of them to be positive, which 
sets an double-sided bound on A^. 

We also studied the relationship between the Yukawa couplings A and k, and we found 
no correlation between them. Even for A = 1.5, the value of n can still vary from 0.3 to 1. 

Next, we take a look at the Higgs boson masses allowed by the constraints. Since a large 
A can enhance the lightest CP-even Higgs boson mass and thus avoid the little hierarchy 
problem, it is interesting to look at the dependence of the Higgs boson masses on the 
parameter A. 

In FigsJS] and H] we show our scan results in A versus rrih plane and A versus iria plane 
with nth being the lightest CP-even Higgs boson mass and nia the lighter CP-odd Higgs 
boson mass. From Figj3]one can learn that the upper bound of mh increases with A, which 
is expected from Eq.([T]), and for A = 1.5 the value of mh can reach 210 GeV. From FigHlone 
can learn that with the increase of A, a super light CP-odd Higgs boson is gradually ruled 
out, and for A > 1 it is bounded in the range 100 GeV ^ ^ 600 GeV. The properties of 
these Higgs bosons can be quite different from those in the MSSM, and their phenomenology 
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FIG. 3: Same as FigiJ], but for A versus the lightest CP-even Higgs boson mass m/j. 
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FIG. 4: Same as Figl2l but for A versus the Ughtest CP-odd Higgs boson mass nia. 
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at the LHC was discussed in 



28| 



Finally, in order to understand the mechanism used to reproduce the correct dark matter 
abundance, we consider the properties in the neutralino sector. In the NMSSM with large 
tan jS, the component of the lightest neutralino is either higgsino- dominant or bino-dominant 
for a light mass below 80 GeV, but for a heavier mass it is bino-dominant. In FigjSlwe show 
our scan results in the plane of m^o versus A. We see that with the increase of A, the upper 
bound on m^o becomes stringent and eventually it is constrained in the range of 50 ~ 100 
GeV. About the next lightest neutralino X2 found that its mass is constrained in the 
range of 100 ~ 160 GeV for A > 1.2. In order to figure out the annihilation mechanism of Xi 
in providing for the dark matter relic density, we compare the masses of X2 ^ with x\ 
in Figini This figure indicates that x^ is significantly heavier than x\- Since in our scan the 
slepton masses are fixed to 200 GeV, also significantly heavier than Xii we conclude that the 
coannihilation of x? with x^ or with a slepton is generally Boltzmann-suppressed and plays 
an unimportant role in accounting for the dark matter relic density. Note that, as shown 
in FigEl there are some samples around the funnel region 2m^o ~ nia and in this case the 
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FIG. 6: Same as FigM but for m^o and nia versus m^o. 
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annihilation of Xi through the s-channel exchange of a hght a becomes dominant 29| 



V. CONCLUSION 



The NMSSM with a large A is an attractive scenario since it can push up the upper 
bound on the SM-like Higgs boson mass to solve the little hierarchy problem. We examined 
the current experimental constraints on this scenario, which include the direct experimental 
bounds, the indirect constraints from precision electroweak measurements, the cosmic dark 
matter relic density, the muon anomalous magnetic moment, as well as the stability of the 
Higgs potential. Our results showed that for a large A the parameter space is severely 
constrained. For example, for a smuon mass of 200 (500) GeV the parameter space with 
A ^ 1.5(0.6) is excluded, and for A = 1 the allowed ranges are 2.5 ~ 4 for tan/3, 400 ~ 800 
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GeV for Ma, 150 ~ 250 GeV for /i, 150 ~ 300 GeV for Ms and ~ 600 GeV for A^. 

Finally, we would like to point out that our conclusion may be qualitatively applicable 
to other NMSSM-like models such as the Minimal Nonminimal Supersymmetric Standard 



Model (MNMSSM) jSOj], which has similar structure with the NMSSM and can be viewed 
as the low energy realization of the Fat Higgs model [7i]. For example, it has been pointed 
out that for any singlet extensions of the MSSM, regardless the form of its superpotential, 
a large A is always accompanied by a small tan/3 fl. This property, as shown m our paper, 
can either limit the smuon mass or limit A if we require the theory to explain the deviation 
of the muon anomalous magnetic momentum. Another example is about the constraint 
from dark matter. In the MNMSSM we expect that the constraint can limit the relevant 
parameters in a more stringent way than in the NMSSM since the neutralino sector in the 
MNMSSM is exactly same as in the NMSSM but with fixed k = [Ml- 
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